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ABSTRACT 

Using Hubble Space Telescope Advanced Camera for Surveys (HST/ACS) and Wide 
Field Camera 3 (WFC3) observations from the Panchromatic Hubble Andromeda 
Treasury (PHAT), we present new period-luminosity relations for Cepheid variables 
in M31. Cepheid from several ground-based studies are identified in the PHAT pho¬ 
tometry to derive new Period-Luminosity and Wesenheit Period-Luminosity relations 
in the NIR and visual filters. We derive a distance modulus to M31 of 24.51 ± 0.08 in 
the IR bands and 24.32 ± 0.09 in the visual bands, including the first PL relations in 
the F475W and F814W filters for M31. Our derived visual and IR distance moduli dis¬ 
agree at slightly more than a 1-cr level. Differences in the Period-Luminosity relations 
between ground-based and HST observations are investigated for a subset of Cepheids. 
We find a significant discrepancy between ground-based and HST Period-Luminosity 
relations with the same Cepheids, suggesting adverse effects from photometric contam¬ 
ination in ground-based Cepheid observations. Additionally, a statistically significant 
radial trend in the PL relation is found which does not appear to be explained by 
metallicity. 


1 INTRODUCTION 

The pulsating class of variable stars known as Cepheids has 
been studied for a number of reasons, the most prominent 
of which is the existence of a relation between their peri¬ 
ods and luminosities. The Period-Luminosity (PL) relation 
for Cepheid variables casts these stars as standard candles, 
making them ideal for determining distances. As a result, 
Cepheids play a vital role in the cosmological distance lad¬ 
der, being observable from the Local Group to distances of 
tens of Mpc, and overlapping with secondary distance indi¬ 
cators. 

More accurate distances from the Cepheid PL Relation 
lead to improved calibration of stellar luminosities, con¬ 
straints on stellar population synthesis models, and mea¬ 
surements of the Hubble Constant (Hq). Improved precision 
of the Period-Luminosity relation has become increasingly 
important in light of the tension between the Cepheid-based 
determination of Hq and that from the Planck mission (Riess 
et al. 2011, Planck Collaboration 2013). 

The possibility of using M31 as a local anchor for the 
Period-Luminosity relation in the future could benefit the 
effort to improve the distance scale. Although the further 
distance degrades accuracy of measurements compared to 
more local distance estimates (he., Milky Way and Large 
Magellanic Cloud), M31 has the opportunity to diversify 
the methods used to determine distances. Compared to the 
Large Magellanic Cloud, M31 is more akin to external galax¬ 
ies used for even further distance measurements. The use of 


M31 as an anchor could potentially benefit distance deter¬ 
minations for the extragalactic distance scale, and perhaps 
assist in addressing the discrepancy between Hq values. 

The Panchromatic Hubble Andromeda Treasury 
(PHAT) survey of M31 provides a unique opportunity to 
image a significant number of Cepheids previously only ob¬ 
served from ground-based telescopes. Obtaining accurate 
Hubble Space Telescope (HST) magnitudes for these vari¬ 
able stars promises to reduce the error in the PL relation for 
M31, leading the way to a more precise value of the distance. 
Previous ground-based surveys have detected thousands of 
Cepheids in and near the M31 disk {e.g., Stanek et al. 1998; 
An et al. 2004; Kodric et al. 2013), but prior to PHAT, only 
small samples of Cepheids in M31 have been observed with 
HST (Macri et al. 2001; Riess et al. 2012). The extensive 
coverage of the PHAT program opens new opportunities to 
improve the distance to M31 via observations of larger num¬ 
bers of Cepheid variables with HST. 

There are three main contributors to the Cepheid dis¬ 
tance uncertainty: metallicity effects, blending and crowd¬ 
ing, and the uncertainty in the distance to the LMC (Riess 
et al. 2009; Freedman & Madore 2010a). The effort to reduce 
the error and dispersion in the PL relation has made large 
gains in the past two decades (Madore & Freedman 1992, 
Macri 2005, Freedman & Madore 2010), though further im¬ 
provements are on the horizon (Gerke et al. 2011). 

There is as yet no clear consensus on the universality 
of the PL relation and its dependence (Bono et al. 2008), 
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or lack thereof (Majaess 2011), on different passbands or on 
the metallicity of the stellar population. If metallicity does 
affect the PL relation, any uncertainties are likely to be min¬ 
imized in the near-infrared (NIR) bands as compared to the 
visual (Madore & Freedman 1991). Additionally, obtaining 
Cepheid PL relations in the NIR reduces the impact of dust 
on magnitude measurements. Both metallicity effects and 
extinction effects have driven the push towards observations 
of Cepheids at longer wavelengths. 

Blending and crowding can contaminate the photome¬ 
try of Cepheids, making true magnitudes more difficult to 
obtain. This problem is significant for ground-based obser¬ 
vations, which may be biased by up to 0.2 magnitudes by 
blending (Mochejska et al. 2000, Vilardell et al. 2007). How¬ 
ever, point-spread function (PSF) magnitudes obtained from 
HST greatly reduce these effects, thereby reducing the error 
in the measurements of Cepheid magnitudes and thus the 
uncertainty in the distance and Hubble constant as well. 
Furthermore, making use of a Wesenheit magnitude (an in¬ 
dex combining magnitude and color) can help to account for 
reddening variations from star to star by taking individual 
Cepheid colors into account (Madore 1982; Opolski 1983; 
Moffett & Barnes 1986; Madore & Freedman 1991; Caputo 
et al. 2000; Leonard et al. 2003; Ngeow & Kanbur 2005; 
Fiorentino et al. 2007; Bono et al. 2008, 2010; Ngeow 2012). 

Lastly, uncertainty in the LMC distance is thought to 
account for about 5% percent of the uncertainty in the cos¬ 
mological distance scale (Freedman & Madore 2010a; Riess 
et al. 2011; Freedman et al. 2012; Inno et al. 2013). By us¬ 
ing NGC 4258 as an anchor for the distance scale instead of 
the LMC, this uncertainty can be reduced to slightly more 
than 3% (Riess et al. 2009b; Riess et al. 2011; Macri et 
al. 2006; Fiorentino et al. 2013). Humphreys et al. (2013) 
have used masers to update the distance to NGC 4258 to 
7.60±0.17±0.15 Mpc, which reduces the error of using the 
LMC as an anchor as well as reducing the error to NGC 
4258 from previous studies (Herrnstein 1999). 

In this work, we present magnitudes of Cepheids in M31 
from the PHAT survey and use PL relations to redetermine 
the distance modulus of M31. Section presents the de¬ 
tails of the data we have analyzed. The methods and the 
analysis of the data to construct PL relations are discussed 
in Section where we decrease the dispersion in the vi¬ 
sual PL relations and increase the sample of Cepheids with 
Near-Infrared (NIR) photometry with respect to Riess et 
al. (2012) and Kodric et al. (2015). In Section we ex¬ 
amine the benefits from the photometric precision from the 
PHAT survey. As we show by examining a subset of our 
Cepheids which were also observed by the DIRECT survey 
(see Section]^ for details), magnitudes obtained via HST in 
the visual bands are superior to recent ground-based sur¬ 
veys (e.g. Stanek et al. 1998; An et al. 2004; Kodric et al. 
2013). In Section we discuss the determination of dis¬ 
tances from the PHAT photometry. The subset of stars from 
DIRECT is used to compare PHAT determined distances to 
ground-based distance estimates. We explore the relation¬ 
ship between metallicity, distance, and radial location for 
the Cepheids in Section 1^ The implications of our results 
are described in Section [Tjwhile the conclusions are summa¬ 
rized in Section Is] 



Figure 1. The locations of 175 Cepheids from the complete 
dataset are shown as circles and the solid line shows the outline 
of the PHAT footprint and coverage of M31. 


2 DATA 

The data used in this paper are obtained from the PHAT 
photometry of M31. The PHAT survey coverage, design, 
and photometry are described at length in Dalcanton et al. 
(2012). The observations utilize the Wide Field Camera 3 
(WFC3) in both the ultraviolet-visible (UVIS) mode and 
the infrared (IR) mode and the Advanced Camera for Sur¬ 
veys (ACS) Wide Field Channel (WFC). They cover about 
one third of the disk of M31 across 6 filters (F275W and 
F336W in the UV with WFC3; FllOW and F160W in the 
IR with WFC3; F475W and F814W in the visual bands 
with ACS). Each HST pointing is two orbits: 1 orbit for 
WFC3/UVIS and 1 orbit for WFC3/IR (with ACS/WEC 
in parallel mode). The coverage of M31 by PHAT is divided 
into 23 “bricks”, each composed of 18 HST fields of view (in a 
6x3 layout, each brick covering a 3x1.5 kpc area). Photom¬ 
etry is performed on each field and brick using DOLPHOT 
(Dolphin 2000) and the subsequent photometry files are fil¬ 
tered to reject low-quality and non-stellar objects, as de¬ 
scribed in Dalcanton et al. (2012) and Williams et al. (2014). 

To obtain the largest sample of Cepheid observations 
possible, we use multiple Cepheid catalogs: PAndromeda 
(Kodric et al. 2013), DIRECT (Kaluzny et al. 1998, Stanek 
et al. 1998, Kaluzny et al. 1999, Stanek et al. 1999, Bonanos 
et al. 2003), and Cepheids already identified in the PHAT 
images (Riess et al. 2012). 

The DIRECT survey imaged M31 from 1996 to 
1997 with the McGraw Hill Telescope at the Michigan- 
Dartmouth-MIT Observatory (MDM) and from 1996 to 
2000 with the 1.2 m telescope at the F. L. Whipple Obser¬ 
vatory (FLWO). The Cepheids in the DIRECT survey were 
observed more than 130 times in the Johnson V-band but 
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less frequently in the Johnson B and Cousins-Kron I filters. 
All stars have V photometry, but not necessarily B or I mag¬ 
nitudes. 6 fields were imaged by the survey, 5 of which over¬ 
lap with the PH AT coverage of M31, including 87 Cepheids. 
Further details about the acquisition and reduction proce¬ 
dures can be found in the DIRECT papers (Kaluzny et al. 

1998, Stanek et al. 1998, Kaluzny et al. 1999, Stanek et al. 

1999, Bonanos et al. 2003). 

The Pan-STARRS survey of M31, also referred to as 
PAndromeda, uses the 1.8 meter Panoramic Survey Tele¬ 
scope and Rapid Response System with the Giga Pixel Cam¬ 
era in Haleakala, Maui, Hawaii (Hodapp et al. 2004, Kaiser 
et al. 2002, Tonry & Onaka 2009). From mid-2010 to late- 
2011, 183 epochs of data were gathered with a half hour of 
observing each night in the rpi and ipi (where PI refers 
to the filter set used on the Pan-STARRS 1 system). De¬ 
tails of the observations as well as reduction procedures can 
be found in Kodric et al. (2013) and Lee et al. (2012). A 
similar study examining these stars in the NIR filters from 
the PHAT images from MAST can be found in Kodric et 
al. (2015), though with a smaller sample size of about 111 
long-period Cepheids. 

Photometry of 67 Cepheids from the first year of PHAT 
data were examined by Riess et al. (2012). The initial iden¬ 
tification of these Cepheids was through the POMME Sur¬ 
vey (Fliri et al. 2012), with an additional Cepheid from the 
DIRECT Survey (Kaluzny et al. 1998, Stanek et al. 1998, 
Kaluzny et al. 1999, Stanek et al. 1999, Bonanos et al. 2003), 
and two more Cepheids from the PAndromeda Survey (Ko¬ 
dric et al. 2013). 

Using all three sources, we produce a dataset containing 
a total of 175 distinct variables with published periods longer 
than 10 days, detailed in Table and whose locations are 
shown in Figure We restrict our sample of Cepheids to 
those with periods longer than 10 days because the presence 
of a break in linearity in the P-L relation at a period of 
approximately 10 days has been well documented (Simon & 
Lee 1981; Tammann & Reindl 2002; Tammann et al. 2002; 
Kanbur & Ngeow 2004; Sandage et al. 2004; Ngeow et al. 
2005; Kodric et al. 2015). The distribution of these periods, 
separately by publication and altogether, is shown in Figure 

El 

To construct the dataset used in this paper, Cepheid lo¬ 
cations were extracted from the PAndromeda and DIRECT 
surveys as well as from Riess et al. (2012). The positions 
were used to locate the photometry in the PHAT database. 
Because it is the largest dataset, positions were taken from 
the PAndromeda survey first, then from Riess et al. (2012), 
ignoring duplicates from PAndromeda, and lastly from the 
DIRECT survey, ignoring duplicates from the other two sur¬ 
veys. We adopt periods from these surveys in the same se¬ 
quence. 

Because the FWHM of the ground-based photometry is 
significantly larger than that of HST, there can be multiple 
sources that are plausible matches based solely on position. 
To identify the most likely Cepheid counterpart in the HST 
data, we first identify sources within 1” of the published 
right ascension and declination. Of these, the final match 
was based on choosing sources with magnitudes within ±1 
magnitude of the source catalog. Objects were rejected if 
they did not show variability beyond the photometric er¬ 
ror in two PHAT epochs. Objects were also thrown out if 



LoglPeriod (days)) Log(Period (days)) 


Figure 2. The original period distributions of each (non cr- 
clipped) sample. Panel a) shows the complete dataset; panel b) 
gives the period distribution for PAndromeda Cepheids; panel c) 
shows the Riess et al. (2012) sample of periods; panel d) gives 
the period distribution of the DIRECT Cepheids. The periods 
are published values from the respective sources seen in Table 


there were multiple objects within 1 magnitude of the ex¬ 
pected, published magnitude(s). For two objects published 
only in the DIRECT dataset, the matching radius was ex¬ 
tended to 1.2” when no match was initially found within 1” 
(the expansion to a larger radii is not surprising given the 
large FWHM of the DIRECT photometry). A comparison of 
the PHAT Cepheids locations to their published locations is 
shown in Figure]^ The mean offsets are approximate 0.11” 
in right ascension and -0.06” in declination. 

Through this method, magnitudes in the PHAT dataset 
were obtained for each Cepheid in each available visual and 
NIR filter. The total sample is composed of 175 Cepheid 
variables with visual and 174 Cepheids with NIR magni¬ 
tudes (see Table one Cepheid is on the very edge of the 
PHAT coverage and falls outside the WFC3/IR footprint). 
The locations of the Cepheids as determined by PHAT are 
given in columns 1 and 2 of Table In columns 3-10 of 
the table, we give the PHAT magnitudes in the visual filters 
(F475W and F814W) and the infrared filters (FllOW and 
F160W), along with their photometric uncertainties from 
photo counts (uncertainties can be significantly larger due 
to crowding; see Dalcanton et al. 2012 ). In column 11, we 
give the published period of each Cepheid as drawn from by 
the reference listed in the final column. 

The magnitudes obtained from the PHAT photometry 
are not time-averaged mean magnitudes, due to insufficient 
temporal coverage. These magnitudes are therefore snap¬ 
shots of the Cepheid at a single phase in its variation. These 
“random” phase magnitudes may lead to deviations from a 
given PL relation due to the amplitudes of the Cepheids’ 
variation and the phase of observation. 

The complete dataset is used for the primary analy¬ 
sis of this paper. However, as we describe in Section]^ we 
separately analyze the DIRECT sample to examine differ¬ 
ences in the PL relations resulting from ground-based and 
space-based photometry. The properties of the DIRECT 
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ARA (arcseconds) 

Figure 3. A comparison of right ascension and declination in the 
PHAT photometry to the published values of the PAndromeda, 
Riess et al. (2012), and DIRECT surveys (i.e.: RAp^fy^j’- 
RAsuruet/). Green circles denote the comparison to PAndromeda, 
blue squares denote the comparison to the DIRECT survey, and 
the red triangles show the comparison to Riess et al. (2012). Sev¬ 
eral Cepheids only found in the DIRECT et al. dataset are outside 
the initial matching radius of 1”; no match was found until the 
matching radius was expanded to 1.2”, allowing a unique match. 
We find a mean offset of 0.11” in right ascension and —0.06” in 
declination. 

sub-sample is seen in Table The DIRECT variable name 
is given in column 1, with the right ascension and declina¬ 
tion as determined by the PHAT data given in columns 2 
and 3. The F475W and F814W magnitudes from the PHAT 
survey are given in columns 4 and 6 with their correspond¬ 
ing photometric errors (columns 5 and 7). Column 8 gives 
the period published by the DIRECT survey. Out of the 85 
Cepheids in the DIRECT sample whose positions overlap 
with the PHAT footprints, only 80 are used due to match¬ 
ing problems, equivalent to a 6% rejection rate. 

2.1 Uncertainty Determination 

We use artificial star tests to determine the true photomet¬ 
ric uncertainties in all 4 relevant filters. For each field and 
camera in the PHAT survey, 10® artificial stars have been 
inserted and reprocessed through the PHAT pipeline (as de¬ 
tailed in Dalcanton et al. 2012 and Williams et al. 2014). Ar¬ 
tificial stars are inserted individually and the photometry is 
re-run in the immediate vicinity of that star and the result¬ 
ing magnitude is recorded. This procedure allow us to fully 
characterize non-trival noise from the photometric measure¬ 
ments, which in particular includes blends and completeness 
From the artificial star tests for the corresponding brick 
and field for each Cepheid, we choose artificial stars with a 
recovered magnitude within 0.5 magnitudes of the Cepheid’s 
observed magnitude and locations within 20 arcseconds of 
the Cepheid’s location. The estimated dispersion and sys- 
tematics of our measurements are then added in quadra¬ 


ture to produce the total photometric uncertainty for each 
Cepheid. 

2.2 Comparison with Riess et al. (2012) 

There are 67 Cepheids in Riess et al. (2012) for which we 
have independently obtained photometry from the PHAT 
photometry pipeline. We compare the FllOW and F160W 
magnitudes for each of these Cepheids in Figure]^ We find a 
median offset of 0.259T0.028 in FllOW and 0.035±0.010 in 
F160W. These offsets are similar to those noted by Kodric 
et al. (2015). Our F160W photometry and that of Kodric 
et al. (2015) differ by an average of 0.016 and our FllOW 
photometry differs by 0.001. Kodric et al. (2015) and the 
PHAT team both utilize PSF photometry and thus achieve 
very similar results. The discrepant F160W point seen in 
Figurej^is the same Cepheid shown to be a misidentification 
in Kodric et al. (2015). 

Figure|^shows the FllOW and F160W color differences 
for the Riess et al. (2012) photometry, the PHAT photom¬ 
etry, and theoretical colors from a grid of isochrones. Our 
set of models is composed of Girardi isochrones over an age 
range of 4 Myr to 1 Gyr and metallicity range of Z=0.0001 to 
Z=0.05, constrained to the canonical instability strip defined 
in Bono et al. (2005). The Riess et al. photometry is shown 
in the top panel and the PHAT photometry in the middle 
panel; the offset between the two observed datasets is clear. 
The offset is due to an ensquared energy fraction correction 
rather than an encircled energy fraction correction, the lat¬ 
ter of which is the correct application for PSF photometry 
(Kodric et al. 2015). Whereas the ensquared energy fraction 
gives the energy over a certain number of pixels, the encir¬ 
cled energy fraction is the light within a certain radii. The 
difference between using the two energy fraction corrections 
is ~0.258, which accounts for the majority of the discrep¬ 
ancies in photometry we see between our photometry and 
Riess et al. (2012). The range and mean of FllOW—F160W 
colors observed in the PHAT photometry mimic those seen 
in the theoretical colors in the bottom panel of Figure 

It is important to note this offset between Riess et al. 
(2012) and our photometry, as it informs the process by 
which we determine a distance with the NIR filters (Section 



3 PERIOD-LUMINOSITY RELATIONS 
FROM PHAT PHOTOMETRY 

We construct optical PL relations using the F475W and 
F814W filters and NIR PL relations using the FllOW and 
F160W filters. The measurements used in these relations 
are “random phase” magnitudes rather than time-averaged 
mean magnitudes. This limitation will introduce scatter into 
the PL relation. Although mean magnitudes lead to less 
scatter in the PL relations, the PHAT data alone does not 
have enough temporal coverage to derive mean magnitudes, 
unlike optimized ground-based Cepheid surveys. 

To determine whether we can accurately recover the 
mean PL relation and avoid biased distance estimates even 
with random phase magnitudes, we did Monte Carlo simu¬ 
lations of 165 light curve random samples to simulate our 
Cepheids using templates from Pejcha & Kochanek (2012). 
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Table 1. Complete Cepheid Sample 


RA 

Dec 

F475W 

^F475W 

F814W 

<7FS14W 

FllOW 

(7F110W 

F160W 

e'Fieow 

Period (days) 

Source 

10.91937 

41.21261 

22.948 

0.052 

20.682 

0.040 

19.982 

0.073 

19.303 

0.100 

10.045 

PAndromeda 

11.27651 

41.69481 

20.764 

0.037 

19.608 

0.053 

19.333 

0.057 

18.883 

0.046 

10.054 

PAndromeda 

11.01628 

41.62777 

21.076 

0.009 

19.796 

0.033 

19.213 

0.081 

18.786 

0.088 

10.101 

PAndromeda 

11.46423 

42.14078 

21.448 

0.006 

19.823 

0.004 

19.173 

0.052 

18.719 

0.063 

10.234 

PAndromeda 

11.02203 

41.23451 

22.039 

0.036 

20.254 

0.025 

19.686 

0.105 

19.134 

0.111 

10.29 

DIRECT 

11.09179 

41.35495 

21.736 

0.014 

19.940 

0.021 

19.490 

0.101 

18.989 

0.114 

10.296 

PAndromeda 

11.18490 

41.92719 

21.163 

0.005 

19.796 

0.017 

19.452 

0.097 

18.895 

0.108 

10.3 

Riess 

11.17153 

41.40722 

22.208 

0.037 

20.239 

0.024 

19.399 

0.100 

18.858 

0.105 

10.35 

DIRECT 

11.34817 

42.03007 

21.452 

0.034 

19.829 

0.017 

19.666 

0.027 

19.076 

0.038 

10.371 

PAndromeda 

11.38078 

41.88077 

21.902 

0.035 

20.033 

0.020 

19.530 

0.044 

19.056 

0.070 

10.43 

Riess 


Table 2. DIRECT Cepheid Sample 


DIRECT ID 

RA 

Dec 

F475W 

(LFA75W 

F814W 

fTpSlAW 

Period (days) 

V5343 

11.02203 

41.23451 

22.039 

0.036 

20.254 

0.025 

10.29 

V8515 

11.09179 

41.35495 

21.736 

0.014 

19.940 

0.021 

10.308 

V13153 

11.17153 

41.40722 

22.208 

0.037 

20.239 

0.024 

10.35 

V2293 

11.13846 

41.62605 

21.120 

0.004 

19.548 

0.006 

10.567 

V6363 

11.36852 

41.65947 

22.249 

0.008 

20.258 

0.007 

10.593 

V410 

11.09179 

41.66419 

21.800 

0.035 

19.999 

0.020 

10.792 

V13042 

11.16980 

41.39944 

21.772 

0.035 

20.141 

0.022 

10.847 

V3773 

10.98764 

41.24036 

23.130 

0.070 

19.874 

0.020 

10.938 

V7381 

11.09263 

41.32135 

20.808 

0.029 

19.586 

0.019 

10.943 

V4733 

11.33250 

41.78882 

22.337 

0.006 

20.241 

0.006 

10.971 



Figure 4. A comparison of the magnitudes from PHAT photom¬ 
etry and Riess et al. (2012) photometry, with FllOW denoted 
as green circles and F160W as blue triangles. The median dif¬ 
ference in the FllOW photometry is 0.259ib0.028 magnitudes, 
a non-negligible offset. The median difference in F160W is only 
0.035ib0.010 magnitudes. 


We then took the mean magnitude (for both B and I) of 
these 165 random samples and compared it to the true mean 
magnitude of the template. This test was repeated 100 times 
and the results are shown in Figure The mean difference 
between the 165 random samples and the true magnitude is 


0.0003 for B and -0.0001 for I; the difference for the B-I color 
is 0.0004. While there may be additional scatter around the 
mean due to sampling a set of Cepheids at random phase, 
there is no apparent bias of the mean in either magnitude 
or color. This suggests that random phase observations do 
not bias our mean distances. 


3.1 Optical Period-Luminosity Relation 

We construct an optical Period-Luminosity relation using 
the random phase magnitudes from PHAT for the Cepheids 
in Table [T] We determine a Wesenheit index for the F475W 
and F814W filters using extinction parameters from the ap¬ 
pendix of Schlafly & Finkbeiner (2011) with an Rv =3.1 
extinction law. The Wesenheit index takes both the bright¬ 
ness and color of each Cepheid into account; this index is 
commonly used in variable studies to correct for star-to-star 
differences in extinction and typically reduces the overall 
scatter in the P-L relation. 


WFi75w,FsiAw ~ FSUW — 0.879(1^47514^ — 1^81414^) (1) 

We use iterative 2.5 -ct clipping with respect to the We- 
senheit(F475W, F814W) relation for the 175 Cepheids with 
F475W and F814W magnitudes, leaving the final sample 
with 163 Cepheids, rejecting 7%. We employ an iteratively 
re-weighted least squares, allowing the slope to float during 
the fitting and clipping process as opposed to clipping with 
respect to a fixed slope. As we sigma-clip with respect to the 
Wesenheit relation rather than the individual bands, we do 
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F110W-F160W 

Figure 5. We compare the FllOW—F160W colors of Riess et 
al. (2012) (top panel), our PHAT photometry (middle panel), 
and theoretical colors (bottom panel). The theoretical colors are 
derived from a grid of Girard! isochrones over an age range of 
4 Myr to 1 Gyr and metallicity range of Z=0.0001 to Z=0.05, 
constrained to the canonical instability strip defined in Bono et 
al. (2005). The mean of the Riess et al. IR colors is 0.32; for the 
PHAT colors the mean is 0.53; the mean derived from isochrones 
is 0.48. The theoretical FllOW—F160W colors derived from the 
grid of isochrones is much more closely matched in mean and in 
range by our PHAT photometry than by the Riess et al. (2012) 
photometry. 



Mean A 


Figure 6. Results of 100 Monte Carlo simulations of 165 ran¬ 
dom light curve random samples compared to the true template 
mean. The top panel shows the mean difference of the 165 random 
sample and the true magnitude for the B filter, the middle panel 
shows the same for the I filter, and the bottom panel shows the 
results for the B—I color. The mean difference of the 165 random 
samples and the true magnitude is 0.0003 for B and -0.0001 for I; 
the difference for the B-I color is 0.0004. While there is scatter, 
there is no apparent bias in the mean. 


not see a bias flattening the slope in the higher-dispersion 
relations. The PL relations are seen in the right panel of 
Figure[7] Table|^gives the slopes and intercepts of the error- 
weighted linear regression relations for F475W, F814W, 
and the Wesenheit magnitudes (WF 475 ty,F 8 i 4 vy, as defined 
in Equation [^. Most notably, the dispersion in the Wesen- 
heit(F475W, F814W) relation is only ~0.17 mag, which is 
the lowest published dispersion for an optical PL relation in 
M31. Moreover, this dispersion is comparable to that found 
in the NIR bands (see Section 3.21. 

However, due to the lack of published PL relations 
in HST native Alters, distance determinations are difficult. 
Thus, we convert the F475W and F814W magnitudes to B 
and I magnitudes (respectively) in order to find distances 
from the optical PL relations. We use the relations from Ta¬ 
ble 2 of Saha et al. (2011) to determine B and I magnitudes. 
We then calculate the Period-Luminosity relationship using 
the Wesenheit magnitude system. We adopt 


Wbi = B- 1.866(5 - I) 


( 2 ) 


from Fouque et al. (2007) with an Rv =3.1 reddening law. 

Note that the color multiplier in the Wesenheit calcu¬ 
lation is a function of the ratio of extinction in the I and 
B filters. However, our B magnitudes are magnitudes from 
the F475W passband converted to B, using semi-empirical 
relations that are functions of stellar temperature and typi¬ 
cally very low extinction. To measure the impact of varying 
extinction on the accuracy of the Wesenheit relation using 
B and I magnitudes transformed from F475W and F814W 
magnitudes, we used Girardi et al. (2008) isochrones, based 
on synthetic spectra, to compare the Wesenheit magnitudes 
at different extinctions (Av=0 and Av = l) with Rv=3.1. 

To begin, we identified points in solar-metallicity 
isochrones (Girardi et al. 2000, Marigo et al. 2008) falling 
in the instability strip in F475W and F814W magnitudes. 
The F475W and F814W magnitudes were converted to B 
and I magnitudes using the Saha et al. (2011) relations. We 
used the B and I magnitudes to compute Wesenheit magni¬ 
tudes in the same way as our observed Cepheid magnitudes 
for Av=0 and Av=l. This process was repeated for We¬ 
senheit magnitudes computed directly from B and I mag¬ 
nitudes from the same isochrone points inside the instabil¬ 
ity strip. We found that our Wesenheit BI magnitudes are 
not extinction-independent, but instead have a dispersion of 
0.06 magnitudes between Av=0 and = In contrast, the 
same experiment done with synthetic B and I magnitudes 
shows a dispersion of 0.01 magnitudes. We fold an error of 
0.06 into our Wesenheit magnitude uncertainty to reflect the 
error associated with making the filter transformation. 

For the 175 Gepheids with visual magnitudes, we use a 
2.5-0- clipping with respect to the Wesenheit relation, leav¬ 
ing the final sample with 163 Gepheids seen in Figure 
rejecting 7%. As with the native filters, we allow the slope 
to float while fitting the PL relations and do not observe a 
flattening of the slope in the noisier relations. Table gives 
the slopes and intercepts of the error-weighted linear regres¬ 
sion relations for 5, 7, and the Wesenheit magnitudes (W bi, 
as defined in Equation]^. The dispersion in the Wesenheit 
BI relation is only 0.19 mag, which is comparable to the 
dispersion in the Wf4,75w,fsi4,w relation and also similar to 
the dispersion in the NIR filters (see Section 3.2 I. 
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Our dispersion is less than other visual PLW (Wesen- 
heit) relations from recent studies (e.g.: >0.34 mag from the 
DIRECT survey, 0.33 mag from Kodric et al. 2013), thereby 
reducing the uncertainty in the eventual distance determi¬ 
nation. Additionally, while the sample size is smaller than 
many ground-based surveys, it is the largest sample size of 
Cepheids observed in M31 with HST in the visual bands. 
The high quality of HST photometry significantly reduces 
the dispersion seen in the PL and PLW relations, attributed 
primarily to the reduction of blending and crowding with 
HST observations. Together, these two improvements lead 
to gains in constraining the uncertainty in the distance mod¬ 
ulus to M31. 

We find a slope of -3.33±0.09 for the Wsr relation, 
which is shallower than the slope found by Fouque et al. 
(2007) with Galactic Cepheids (-3.600T0.079) and their 
slope with LMC Cepheids (-3.454T0.011). However, Fouque 
et al. (2007) includes a period range from logj^Q(P) « 0.6 to 
1.8 in determining a relation, which may contribute to some 
discrepancy, as including shorter period Cepheids tends to 
steepen the PL slope. 


3.2 Infrared Period-Luminosity Relation 

In Figure we show the Period-Luminosity relation from 
the random phase magnitudes for the Cepheids in FllOW 
and F160W WFC3 filters, and in Wesenheit magnitudes, 
as calculated as in Riess et al. (2012) with an Rv =3.1 
reddening law: 


W7fl = F160W-1.54(F110W-F160W) (3) 

As with the optical bands, we fit linear, error-weighted 
relations to the data. The resulting slopes, intercepts, and 
dispersions are given in Table Out of the 175 Cepheids 
with IR magnitudes, 160 are left after 2.5-(7 clipping with re¬ 
spect to the Wesenheit relation (allowing the slope to float, 
producing a 9% rejection). Although we present the NIR 
Period-Luminosity relations in Table we do not directly 
use these fits to determine a distance in Section [5.21 As ex¬ 
pected, the individual NIR bands give a much smaller scatter 
in the Period-Luminosity relation, due to lower extinction 
and smaller overall amplitudes of the Cepheid light curves. 
However, although the amplitudes in the NIR are generally 
less than those in the visual, they are not completely negli¬ 
gible (Monson & Pierce 2011, Testa et al. 2007, Persson & 
Madore 2004, Welch et al. 1984). 

The dispersion in the Wesenheit relation is similar to 
previous studies and comparable to that of the visual Wesen¬ 
heit dispersion. Our results show slightly higher dispersion 
in FllOW (0.246 mag) than found by Riess et al. (2012, 0.20 
mag), but very similar to that of Kodric et al. (2015, 0.243). 
It is possible the discrepancy in the dispersions between our 
study and Riess et al. (2012) is related to the differences in 
photometry. Similar dispersions in the F160W PL are found 
in all three studies: 0.187 mag in this paper; 0.17 mag in 
Riess et al. 2012; 0.178 in Kodric et al. (2015). We find a 
dispersion of 0.173 mag in the PLW relation, between 0.22 
from Riess et al. (2012) and 0.147 Kodric et al. (2015), al¬ 
though we note that Kodric et al. use a different dehnition 
of the Wesenheit magnitude. Save for FllOW, our disper¬ 
sions are also consistent with the expected random phase 



Figure 8. NIR Period-Luminosity relations for 160 of the 175 
Cepheids with magnitudes in FllOW and F160W and Wj_r. The 
F160W magnitudes are shifted brighter by 2 magnitudes and the 
Wesenheit magnitudes are shifted brighter by 4 magnitudes to 
make the three relations clearly visible. Weighted least-squares 
linear fits to each relation after 2.5-cr cuts are plotted on top 
of the data. Uncertainties determined by artificial star tests are 
included, though for some points they are too small to see clearly. 


dispersions as noted by Riess et al. (2012) from the study of 
Persson et al. (2004) of LMC Cepheids. 

Generally, the IR bands provide an opportunity to im¬ 
prove upon the dispersion of the PL relation typically seen in 
the visual bands. However, we do not necessarily see a signif¬ 
icant difference between the dispersions of PLW relations in 
the visual and in the IR, although the individual bands show 
immense improvement compared to the visual as the wave¬ 
length gets longer, as expected. By increasing the sample 
size of long-period Cepheids in the IR bands over Riess et al. 
(2012, 68 Cepheids) and Kodric et al. (2015, 110 Cepheids), 
we have further reduced the statistical uncertainty in the 
distance measurement. 

Kodric: FllOW: -2.497±0.209 F160W: -2.779±0.171 

Comparing the linear fits of the F160W filter from Ta¬ 
ble 3.2 to Kodric et al. (2015), we find a similar slope (- 
2.96±0.09 from this paper and -2.779±0.171) and intercept 
(19.02±0.02 from this paper and 18.960±0.028 from Kodric 
et al. 2015). The FllOW slope from Kodric et al. (2015) of 
-2.497T0.209 is significantly flatter than the slope we find 
of-2.92T0.il. The intercepts are reasonably similar within 
1.5-cr (19.58T0.03 from this paper and 19.476T0.037 from 
Kodric et al. 2015). 

We find that our IR Wesenheit slope of -3.38T0.08 is 
consistent with Riess et al. (2012), who found a slope of 
—3.43T0.17 and with Persson et al. (2004) for the LMC, who 
found a slope of -3.38T0.09. The period range of Persson et 
al. (2004) (extending to logjo(P)~2) is comparable to that of 
Riess et al. (2012) and our own complete sample in the NIR 
(extending to logj^Q(P)«1.9, see Figuresand [^. Our slope 
is steeper than that of Kodric et al. (2015) at -3.172T0.117; 
however, they calculate the Wesenheit magnitude slightly 
differently so the comparison is incomplete. 
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Table 3. Optical Period-Luminosity Relations 


Magnitude 

Slope 

(mag/ Log P) 

Intercept (P=10d) 

(mag) 

Dispersion to Fit 

(mag) 

B 

-1.44±0.55 

22.48±0.14 

0.922 

I 

-2.40±0.24 

20.15±0.06 

0.409 

Wesenheit (BI) 

-3.39±0.09 

18.29±0.02 

0.190 

F475W 

-1.62±0.5 

23.72±0.13 

0.833 

F814W 

-2.41±0.24 

22.56±0.06 

0.403 

Wesenheit (F475W,F814W) 

-3.27±0.08 

21.83±0.02 

0.171 



Figure 7. Left panel: optical Period Luminosity relations for 163 of the 175 Cepheids with transformed B and I magnitudes and Ws/. 
The I magnitudes are shifted brighter by 2 magnitudes and the Wesenheit magnitudes are shifted brighter by 3 magnitudes to make the 
three relations clearly visible. Right panel: optical Period Luminosity relations for 163 of the 175 Cepheids with magnitudes in F475W 
and F814W and Fi75W,F6l4W■ The F814W magnitudes are shifted brighter by 2 magnitudes and the Wesenheit magnitudes are shifted 
brighter by 4 magnitudes to make the three relations visible. In both panels, weighted least-squares linear fits to each relation after 2.5 -(t 
cuts are plotted on top of the data. Uncertainties determined by artificial star tests are included, but often too small to see. 


4 GROUND VS. HST PL RELATIONS 

We follow a similar procedure as Section [3.1| to analyze the 
subset of DIRECT Cepheids in the PHAT photometric cat¬ 
alog. From this, a comparison of ground-based and HST 
photometry can be made. We use the PHAT photometry 
and published DIRECT periods and magnitudes to analyze 
the differences between space-based random phase magni¬ 
tudes and ground-based mean magnitudes. For each of the 
80 Cepheids in Tablewe proceed as follows. First, we note 
that, because the DIRECT survey is mainly in the V and I 
(only about half the Cepheids have B magnitudes) bands, 
the U-band and 7-band magnitudes must be inferred from 
the HST data to make a comparison of the PHAT survey to 
the DIRECT survey. 

To transform F814W magnitudes to I magnitudes, we 
use the Saha et al. (2011) relations. We determine a color 
transformation from the F475W and F814W filters to the 
U-band Johnson-Cousin filter through the use of Girardi 
isochrones at a variety of ages (1 Gyr through 10 Gyr) and 
extinctions (Ay from 0 to 2 mags). The transformation is 
restricted to the F475W-F814W color range of our DIRECT 
sample (F475W-F814W from 0.9 to 3.2). The conversion we 
determine is: 


V - F814W ^ 0.104 + 0.541(F475tT - F814W) 
-b0.032(E475W - F8l4Wf 

with an dispersion of 0.03 mags in the transformation. The 
PHAT transformed V and I magnitudes are used to con¬ 
struct PL relations and compared to the PL relations using 
the published V and I magnitudes from the DIRECT sur¬ 
vey. 

We calculate Wesenheit magnitudes for each Cepheid 
via the following relation from Fouque et al. (2007), where 
Wy is the Wesenheit magnitude, V is as in Equation]^ and 
I is transformed via Saha et al. (2011) from F814W: 

Wvi = V-2.55{V - I) (5) 

The PHAT PL relations for V, I, and Wesenheit mag¬ 
nitudes are shown in the right panel of Figure For com¬ 
parison, PL relations from the DIRECT ground-based mean 
magnitudes are shown in the left panel of Figure]^ We use 
2.5 (T-clipping with respect to the Wesenheit magnitudes and 
allow the slope to float. Out of the 80 Cepheids in the sam¬ 
ple, this leaves 77 Cepheids in the DIRECT photometry (a 
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Table 4. Infrared Period-Luminosity Relations 


Magnitude 

Slope 

(mag/ Log P) 

Intercept (P=10d) 

(mag) 

Dispersion to Fit 

(mag) 

FllOW 

-2.92±0.11 

19.58±0.03 

0.246 

F160W 

-2.96±0.09 

19.02±0.02 

0.187 

Wesenheit 

-3.38±0.08 

18.45±0.02 

0.173 


4% rejection rate) and 75 Cepheids in the PHAT photome¬ 
try (a 6% rejection rate). Along with their dispersions, the 
linear, error-weighted fits to the relations in V, I, and We- 
senheit magnitudes from Figure]^ are given in Table 

As illustrated in Figure the ground-based DIRECT 
Period-Luminosity relations have much shallower slopes and 
brighter intercepts compared to the HST PL relations. A 
shallower slope leads to a longer distance modulus; however, 
a brighter intercept leads to a shorter distance modulus. 
These two effects may balance out in distance determina¬ 
tions, but it is evident that ground-based observations are 
susceptible to heavily biased slopes and intercepts and cau¬ 
tion should be taken when drawing conclusions from ground- 
based PL relations. This behavior suggests that blending 
and crowding may be the most significant limitation for 
ground-based Cepheid studies; these issues are improved by 
the spatial resolution of ACS/WFC. We note, however, that 
the level of crowding in ground-based surveys of M31 is com¬ 
parable to that expected for HST surveys of more distance 
galaxies. Therefore, we expect many HST Cepheid studies 
in more distant galaxies may be affected by the same bi¬ 
ases affecting the DIRECT observations of M31 (Stanek & 
Udalski 1999, Mochejska et al. 2000, Chavez et al. 2012). 

As Figure shows, the difference between the slopes 
and the intercepts between the ground-based PL relations 
and the PHAT PL relations grows as we move from the 
V band to the I band. The difference is most stark when 
comparing the Wesenheit slopes and intercepts. This im¬ 
plies that the ground-based bias signihcantly affects both 
the slope and the intercept of the PL and PLW relations, 
thereby jeopardizing accurate distance estimates. The exis¬ 
tence of this bias stresses the need for more HST observa¬ 
tions of Cepheids, even if they are only at random phases. 

The slope of the visual Wesenheit VI relation from the 
PHAT data is in agreement with previous studies. We find 
a slope of -3.30T0.12, comparable to that found by Bene¬ 
dict et al. (2007) and Fouque et al. (2007) for Galactic 
Cepheids (-3.31±0.17 and -3.377±0.023) and that of Udal¬ 
ski et al. (1999) for the LMC (-3.28±0.14), although the 
period ranges differ. The periods in Benedict et al. (2007) 
range from very low period (logjQ(P)=0.4) to longer periods 
(logjg(P) = 1.6), but the data points are extremely sparse 
between logj^Q(P) = l and their longest period Cepheid at 
logj^g(P) = 1.6. The sample in Fouque et al. (2007) is bet¬ 
ter populated in period, with a range from approximately 
logj^g(P)=0.6 to logjQ(P)=1.7. For the Cepheids in the com¬ 
plete dataset, the period coverage is out to log]^Q(P)=1.9, a 
bit farther than Fouque et al. (2007). However, we do not 
include the short period (<10 days) Cepheids in our deter¬ 
mination of the PL relations. 


5 DISTANCE MODULI 

5.1 Distances: Visual 

To determine a distance modulus for the complete dataset in 
the visual bands, we use the Wesenheit BI relation as defined 
by Fouque et al. (2007), calibrated to Galactic Cepheids. 
The distance modulus in these filters is defined by: 

Mw{bi) = -3.600(±0.079) logig(P) - 2.401(±0.023) (6) 
and 

Um31bi = 'mw(Bi) — Mw(bi) (7) 

where M.w(bi) is the magnitude of the Wesenheit BI rela¬ 
tion from Fouque et al. (2007) evaluated at logj^g(P) of 1.26, 
the mean period of our sample, and is the magni¬ 

tude of the Wesenheit BI from our observations, also evalu¬ 
ated at logig(P)=1.26. The result is a distance modulus of 
24.34±0.10. 

We also determine a distance modulus using the LMC- 
based relation from Fouque et al. (2007), using the same 
formalism as Equations and except with a slope of - 
3.454T0.011 and an intercept of 15.928T0.003 (adopting a 
distance to the LMC of 18.486T0.065 as in Riess et al. 2011). 
The result is a distance modulus of 24.31±0.07. We aver¬ 
age our two visual band distance determinations to obtain 
a modulus of 24.32 ± 0.09. 

5.2 Distances: IR 

To obtain a distance modulus to M31 for the NIR Cepheids, 
we tie the distance modulus of M31 to that of NGC 4258, 
as defined in Riess et al. (2009b) and Riess et al. (2011) 
(revisited by Riess et al. 2012 and Fiorentino et al. 2013). 
From Equation 7 of Riess et al. (2009b), we solve for the 
distance modulus to M31 as 

PM31 = M4258 — «P4258 + THw “ bw logio(P — 1) 

-fZwAlogio[0/H] 

where /44258 is the distance to the water maser galaxy NGC 
4258 (7.6 ± 0.17 ± 0.15 Mpc from Humphreys et al. 2013), 
zpi 258 is the intercept of the PL relation for Cepheids in 
NGC 4258, mw are the observed NIR Wesenheit magni¬ 
tudes of the Cepheids, hw is the slope of the global Wesen¬ 
heit NIR PL relation, and P are the periods (>10 days). The 
term ZwA logj^Qp/H] is the metallicity dependence term; 
A logj^oP/H] are the metallicities of the Cepheids relative 
to the LMC and Zw is the global metallicity slope. To ob¬ 
tain the values of zp 425 s, bw, and Zw we use a simultaneous 
linear fit to the Cepheids in both M31 and NGC 4258. The 
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Table 5. DIRECT Ground-Based Period-Luminosity Relations 



Magnitude 

Slope 

(mag/ Log P) 

Intercept (P=10d) 

(mag) 

Dispersion to Fit 

(mag) 

DIRECT Ground-Based 

V 

-0.67±0.33 

20.89±0.09 

0.508 


1 

-1.32±0.21 

19.76±0.06 

0.339 


Wesenheit 

-2.33±0.30 

18.01±0.08 

0.440 

PHAT Random Phase 

V 

-1.22±0.44 

21.23±0.11 

0.541 


1 

-2.03±0.24 

20.13±0.06 

0.327 


Wesenheit 

-3.30±0.12 

18.31±0.03 

0.187 



Figure 9. Optical PL relations for ground-based DIRECT (left), for random phase HST (right) DIRECT Cepheids with V, I (transformed 
from F475W and F814W), and Wy/ magnitudes. The I magnitudes are shifted brighter by 2 magnitudes and the Wesenheit magnitudes 
are shifted brighter by 4 magnitudes to make the three relations clearly visible in each panel. The left panel includes the published 
DIRECT ground-based mean magnitudes for each Cepheid, with 77 Cepheids and error bars as given in the published DIRECT data. 
The right panel shows the same Cepheids’ random phase magnitudes determined from the PHAT photometry, with 75 Cepheids and 
error bars as determined from artificial star tests. Weighted least-squares linear fits to each relation after 2.5-cr cuts are plotted on top 
of the data in both panels. 


simultaneous linear fit allows a reduction in the final dis¬ 
tance uncertainty compared to using the NIR PL relations 
alone. 

First, we should note that same filters were not used 
to observe Cepheid colors in NGC 4258 (V-I) as the bands 
observed in the PHAT survey (F160W-F110W). The We¬ 
senheit magnitudes in Riess et al. (2009 and 2011) were cal¬ 
culated using W/i{=F160W-1.54*(V-I). To account for the 
difference in defining the Wesenheit magnitude, a correction 
has been applied to the (F110W-F160W) based extinction 
correction to give them the same mean color in (V-I) as M31 
Cepheids. As in Riess et al. (2012), we use the following to 
determine a correction (X) to the Wesenheit magnitude cal¬ 
culation of NGC 4258: 

0.504(V - 7) = 1.54(7^1101^ - FimW - X) (9) 

With a mean (V-I) color of 1.23 (from the DIRECT 
Cepheids) and a mean (F110W-F160W) color of 0.56 from 
the PHAT photometry, we obtain a value of X=0.156. Be¬ 
cause of the difference in FllOW photometry between PHAT 
and Riess et al. (2012), we use a different X than Riess et 
al. (2012) to place our Wesenheit magnitudes on the same 


scale as that of the NGC 4258 photometry. Errors in their 
NGC 4258 photometry are incorporated into the analysis 
and error budget. 

To solve for the distance, we also must determine the 
values of A log jq [0/77] for the Cepheids to include the possi¬ 
bility that differing metallicities of Cepheids could introduce 
scatter into the P-L relation. To examine possible metallic- 
ity effects, we use the relation from Zaritsky et al. (1994) 
for M31 to determine log]^Q[0/H] for each Cepheid from its 
radial location in the disk and compare it to the solar value 
of 8.9 (as in Riess et al. 2012, Riess et al. 2011, and Riess et 
al. 2009). 

We adopt a center for M31 of RA = 0^ 42"* 44.31® and 
Dec = 41° 16’ 9.4” (Cotton et al. 1999), a disk scale length of 
Po=77.44 (Zaritsky et al. 1994), and position angle (/=37.715 
(Hand et al. 1981, Baade & Arp 1964). These values are 
used to calculate the 12-|-logjQ[0/H] metallicity term from 
the relation determined by Zaritsky et al. (1994), as seen in 
Equation |10[ 


12 + logio([0/77]) = 
9.03(±0.09) - 0.28(±0.10)(p/po - 0.4) 


© 2014 RAS, MNRAS 000, [^?? 

















Panchromatic Hubble Andromeda Treasury XIII: The Cepheid Period-Luminosity Relation in MSI 


11 


DIRECT Ground-Based data 
PHAT {DIRECT Cepheids) 
PHAT (Complete dataset) 



Wesenheit 



Figure 10. Left panel: a comparison of intercept values across several filters; circles denote intercepts from the DIRECT ground-based 
PL relation, triangles indicate the intercepts from the DIRECT sample as observed by PHAT, and squares indicate the intercepts from 
the complete PHAT sample (/ band only). The left section shows intercepts from the I/-band PL relation, the middle from the 7-band, 
and the right section from the W\// PL relation. Error bars are also plot for each intercept. Right panel: a comparison of slope values 
across several filters; circles denote slopes from the DIRECT ground-based PL relation, triangles indicate the slopes from the DIRECT 
sample as observed by PHAT, and squares indicate the slopes from the complete PHAT sample (7 band only). The left section shows 
slopes from the V-band PL relation, the middle from the 7-band, and the right section from the Wyi PL relation. Error bars are also 
plot for each slope. 


The deproiected radius, p, is calculated using Equations 
[T3 through!^ 


X = (S — 5o)cos{(f)) + (a — ao)sin{4>)cos{5o) (11) 


(5 — 5o)sin{(f>) — (a — ao)cos{(l>)cos{5o) 
cos{i) 


( 12 ) 


P = 


po 


(13) 


The resulting 12+logio[0/H] values are compared to 
the IR Wesenheit magnitude in Figure IE The slope of a 
linear fit relating the metallicity of the Cepheids to their 
magnitudes, is 0.22±0.60 mag dex“^. A two-sided student 
t-test with a value of t=0.37 and 159 degrees of freedom 
informs that this is not a statistically significant relation 
at the 0.05 significance level. Our slope of 0.22±0.60 mag 
dex“^ is different than Riess et al. (2012), who obtained a 
slope of -0.65±0.73, but with large uncertainty. 

We see a range of ~0.3 dex in 12-|-log,^|3[0/H] from ap¬ 
proximately 8.82 to 9.12, similar to Riess et al. (2012), who 
find a range of 8.87 to 9.05. This range is similar to the 
~0.4 dex that we expect from the -0.018 dex/kpc gradi¬ 
ent of Zaritsky et al. (1994) over the range of the disk. The 
Cepheids in M31 do not have a clear magnitude-metallicity 
dependence, as has been found in previous studies (Riess et 
al. 2012, Freedman & Madore 1990). However, the values of 
12-|-log,^Q[0/H] for M31 Cepheids will be used in conjunction 
with those in NGC 4258 to find a global fit. 

Using multiple error-weighted linear regression, we 
simultaneously fit a period-luminosity relation to the 



Figure 11. The 12-|-log,^Q[0/H] vs Wm magnitude relation for 
160 Cepheids in M31. The black points indicate the individual 
Cepheids’ metallicities determined by the relation from Zaritsky 
et al. (1994), Equation |10[ with propagated errors. The solid line 
shows the linear, error-weighted fit to the Cepheids, with a slope 
of 0.22±0.60 mag/dex, which is not statistically significant at the 
0.05 level with a two-sided student t-test value of 0.37 and 159 
degrees of freedom. We see a range of ~0.3 dex in 12+logjQ[0/H] 
from approximately 8.82 to 9.12. 


Cepheids in M31 and those in NGC 4258 to determine val¬ 
ues for zpi 25 s, bw, and Zw- These values are then used to 
determine the distance modulus to M31 as in equation]^ We 
determine a hw slope of-3.24±0.10, 2 P 4258 of 23.26±0.06, 
and Zw of 0.16±0.26. These are comparable to the values 
determined by the global fit in Riess et al. (2011, see Table 
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Table 6. Global Fit Values 


Variable (Eq 


8i 


This Paper 


Riess et al. 2011 


1^4258 

7.6±0.17±0.15 

7.2±0.2±0.3 

ZP4258 

23.26±0.06 

26.32±0.03 (23.10 at P: 

hw 

-3.24±0.10 

-3.21±0.03 

Zw 

0.16±0.26 

-0.10±0.09 



Figure 12. Individual IR distances (described in Section |5.2[ | 
with propagated errors are plotted against period. No correlation 
between the distance modulus and the period of variability is 
seen. The dashed line indicates the mean distance. 


[^. Note that while Riess et al. (2011) use all SN hosts for a 
global fit, we use NGC 4258 with M31 only. With the values 
in Table]^ we obtain a distance of 24.51 ± 0.08 to M31. 

The derived distance moduli for the complete dataset in 
the visual (with MWG and LMG anchors) and in the IR are 
given in Table There is an offset between the visual and 
the IR distances. We discuss this difference in Section 0 be¬ 
low. We plot the individually determined distances derived 
from the IR photometry against period in Figure there 
does not to appear to be a trend with period. Published 
period uncertainties are not available for the DIRECT data 
and the POMME Cepheids from Riess et al. (2012), so we 
cannot fully incorporate period uncertainties into our anal¬ 
ysis. However, if we adopt a reasonably conservative 0.1 day 
period error, the uncertainty in the NIR and visual distances 
is affected by less than 0.01. 


5.3 Ground vs. HST Distances 

For the DIRECT Cepheid dataset, we use the Wesenheit 
VI relation as defined by Fouque et al. (2007), calibrated 
to Galactic Cepheids. We obtain distance moduli for the 
random phase PHAT Cepheid PL Wesenheit relations and 
the ground-based DIRECT PL Wesenheit relations in the 
V and / bands by comparing the magnitude of our PL re¬ 
lation to the calibrated absolute magnitude PL relation at 
the mean period of our sample (P=1.26). The absolute We¬ 
senheit magnitude is given by Fouque et al. (2007) as: 


Mw{vi) = -3.477(±0.074) logio(P) “ 2.414(±0.022) (14) 
The distance modulus is then defined as: 


— ’’nwiyi) — Myy(vi) (15) 

where Mw(yr) is the magnitude of the Wesenheit VI relation 
from Fouque evaluated at log]^Q(P) of 1.26 (the mean period 
of our sample) and mw(vi) is the magnitude of the We¬ 
senheit VI relation from our observations, also evaluated at 
logio(P)=l-26. We also determine a distance modulus using 
the relation form Udalski et al. (1999), using the same for¬ 
malism as Equations |14| and |15| with a slope of-3.320±0.014 
and an intercept of 15.880±0.010 (adopting an LMC dis¬ 
tance modulus of 18.486±0.065 as in Riess et al. 2011). 

Our W(VI) DIRECT sample extends to a maximum pe¬ 
riod of about logiQ(P)=1.75, comparable to that of Fouque 
et al. (2007). Udalski et al. (1999) have period coverage from 
logio(P)=0.4 to logjQ(P)=1.5, where the LMC Cepheids be¬ 
gin to saturate their detector. In both cases, the authors in¬ 
clude the shorter period Cepheids in determining a PL rela¬ 
tion; however, we restrict our sample to only those Cepheids 
with periods greater than 10 days. 

Distance estimates for the DIRECT subset from both 
the ground-based DIRECT photometry and the PHAT pho¬ 
tometry are included in Table[^ Distances are given for both 
a Milky Way and an LMC anchor. 

The ground-based distance moduli appear to be shorter 
than their HST counterparts by ~0.03. While this differ¬ 
ence is well within the error, it could also be due to less 
accurate photometry on the ground. Blending and crowding 
could plausibly cause systematically shorter distances from 
ground-based photometry. We discuss this issue further in 
section [T] 


6 METALLICITY AND THE 

DISTANCE-RADIUS RELATION 

A slight trend of Cepheid distance modulus with radial loca¬ 
tion from the center of the galaxy has been seen in MlOl and 
M81, among others (see Freedman & Madore 1990, Kenni- 
cutt et al. 1998, Macri et al. 2006, McCommas et al. 2009, 
Gerke et al. 2011). The discrepancy in distance between 
Cepheids in the inner and outer regions of these galaxies 
is often attributed, at least in part, to metallicity effects on 
the Period-Luminosity Relation, although the cause is still 
debated (Gould 1994; Vilardell et al. 2007, Gerke et al. 2011, 
Majaess et al. 2011, Kudritzki et al. 2012). 

We use our complete set of Cepheids from Tablej^to ex¬ 
plore the relationship between distance and radius in M31, 
determining the radius of each Cepheid by using the right 
ascension and declination to determine its de-projected po¬ 
sition in the galaxy, p, as seen in Equation |13| We convert p 
to kpc using the mean IR distance modulus determined in 
Section The individual distances in the IR bands (pm) 
for each Cepheid are compared to the de-projected radii in 
Figure [T^ A linear fit yields a gradient of 0.007±0.003. The 
two-sided student t-test value for this relation is 2.71, mak¬ 
ing it significant beyond the 99.5% level with 159 degrees of 
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Table 7. Distance Moduli 


Dataset Sample 

Filter 

Anchor 

u 

Complete 

PHAT Visual (B, I) 

Galactic Cepheids 
LMC 

24.34T010 

24.31T0.07 


PHAT NIR (FllOW, F160W) 

NGC 4258 

24.51 ± 0.08 

DIRECT 

Ground-based (V, I) 

Galactic Cepheids 
LMC 

24.21T0.11 

24.21T0.09 


PHAT Visual (V, I) 

Galactic Cepheids 
LMC 

24.24T0.10 

24.24T0.08 


freedom. We apply Chauvenet’s criterion to the data to de¬ 
termine the sensitivity of individual points on the relation. 
We find 6 outlying points with a significance level of 0.05. 
While this shallows the slope slightly, the slope remains sta¬ 
tistically significant. This suggests that there is a statistical 
difference in the distance moduli of individual Cepheids in 
the inner and outer portions of M31. 

We also examine the radial magnitude trend in M31 by 
adjusting each Cepheid’s Wesenheit magnitude based on its 
period, thereby correcting for the PL relations we determine 
in Sections [3^ (this is comparable to examining the residual 
of the PL relation). The result of this exercise is plotted in 
Figure We see that the Cepheids tend to be brighter in 
the inner region of the disk and fainter in the outer parts 
of M31. There is a correlation between magnitude and ra¬ 
dial location with a two-sided student t-test at the 95% sig¬ 
nificance level. We also apply Chauvenet’s criterion to the 
data and find 7 outlying points with a significance level of 
0.05; however, the slope is not significantly affected. The ra¬ 
dial trend in magnitude is causing the change in distance 
modulus; however, the cause of the inwards brightening of 
Cepheids is unclear. 

As previously discussed, we use the relation derived by 
Zaritsky et al. (1994) to calculate the 12-|-logjQ[0/H] val¬ 
ues for each Cepheid as a proxy for the metallicity effects in 
M31. Zaritsky et al. (1994) found a radial metallicity gradi¬ 
ent of-0.018±0.006 from their observations of HIT regions in 
M31, which is quite shallow compared to many of the galax¬ 
ies in that study. This is the gradient employed by this study 
and by Riess et al. (2009, 2011, 2012) in defining their dis¬ 
tance scale. Sanders et al. (2012) find a metallicity gradient 
of -0.0195±0.005 in HII regions comparable to that found 
by Zaritsky; however, they do not find the trend echoed in 
their observations of planetary novae. Several studies have 
examined the globular clusters in M31 to search for a metal¬ 
licity gradient. Barmby et al. (2000) looked at both spectro¬ 
scopic and photometric metallicities of globular clusters, and 
while overall found no trend, observed a slight gradient of 
-0.023±0.01 when restricting their sample to clusters with 
spectroscopic metallicities. Huxor et al. (2011) use CMD 
metallicities but assert the lack of a metallicity gradient, 
suggesting that the gradient was driven primarily by metal- 
poor clusters. Freedman & Madore (1990) use Cepheids in 
3 different fields of M31 and, despite there being a slight 
difference in distance moduli between the inner and outer 
fields, conclude that their data is consistent with there be¬ 
ing no dependence of the PL zeropoint on metallicity. How¬ 
ever, Gould (1994) re-analyzed the same BVRI data used 
by Freedman & Madore (1990) to show the M31 distance 


should be corrected for metallicity to avoid introducing ad¬ 
ditional errors or systematics. Lee et al. (2013) use 17 beat 
Cepheids to trace the metallicity in M31 and find a gradi¬ 
ent of -0.008±0.004, shallower than that of Zaritsky et al. 
(1994) and Sanders et al. (2012), but similar to the PNe gra¬ 
dient from Kwitter et al. (2012) of-0.011 ±0.004). While the 
metallicity gradient in M31 may not be particularly steep, 
there is a wealth of data supporting its existence and direc¬ 
tion. 

Previous empirical studies have suggested metallicity 
as the cause of the distance-radius relationship. However, 
there is disagreement between empirical studies and pulsa¬ 
tion theory on how strong the metallicity effect is, as well 
as whether a higher metal content leads to a brighter or 
fainter Cepheid (Freedman & Madore 1990; Kennicutt et 
al. 1998; Baraffe & Alibert 2001; Maori et al. 2006; Caputo 
2008; Bono et al. 2008, 2010; Romaniello et al. 2008, and ref¬ 
erences therein). Most theoretical approaches predict that 
metal-rich Cepheids should be fainter (Baraffe & Alibert 
2001; Caputo 2008). Conversely, most observational studies, 
including indirect and direct metallicity measurements, sug¬ 
gest that the more metal-rich Cepheids are brighter (Kov- 
tyukh et al. 2005b; Macri et al. 2006; Bono et al. 2008), 
although Romaniello et al. (2008) presented spectroscopic 
evidence based on observations of Cepheid iron lines that 
disagrees with most empirical studies. An alternative vi¬ 
able explanation to metallicity effects could be crowding and 
blending, causing inner Cepheids to appear brighter (leading 
to shorter distances) compared to Cepheids in outer fields. 

Our results show a statistically significant distance- 
radius relationship in the IR bands, but without a significant 
metallicity component, suggesting that there may be an¬ 
other factor at least partially responsible for the difference. 
Majaess et al. (2011) argue that crowding is the cause of 
brighter Cepheids in the inner regions. In comparing ground- 
based and HST, Mochejska et al. (2000) and Chavez et al. 
(2012) both show that blending can be a significant effect, 
especially as seeing increases in the ground-based observa¬ 
tions. Bono et al. (2008) that blending will be a stronger 
effect in the central region of a galaxy and thus lead to a 
decrease in distances. 

Based on our results, crowding and blending could both 
be likely candidates for the cause of the brightening of in¬ 
ner Cepheids in the case of M31. We don’t see a trend with 
magnitude and metallicity, but we do see a trend of mag¬ 
nitude (and distance) with radius. When we examine the 
photometric uncertainties from artificial star tests, we see 
the magnitude error and crowding parameter increase to¬ 
wards the center of M31. The average trend in uncertainty 
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Figure 13. The Distance-Radius Relation for the 160 Cepheids in 
M31 that were used in the final PL relations in Section |3.2| using 
the complete dataset. The black points indicate the individual 
Cepheids’ IR distances with propagated errors. The horizontal 
dotted line shows the mean distance determination from Section 
and the solid line shows the linear fit to the data. The linear 
fit has a slope of 0.007 it 0.003, and a two-sided student t-test 
value of 2.71 (with 159 degrees of freedom), making it significant 
at the 99.5% level. 



Figure 14. The relationship between W/ij magnitude and de- 
projected radius (using the distance modulus determined for the 
IR, /rjij=24.51 it 0.08), where the magnitudes herein are adjusted 
to take into account the periods of the Cepheids, as described in 
the corresponding text, and include propagated errors. The linear 
fit has a slope of 0.005 it 0.002. With a two-sided student t-test 
value of 1.89 and 159 degrees of freedom, the slope is significant 
beyond the 95% level. 


(~0.002 magnitudes/kpc) is sufficient to account for the ob¬ 
served gradient in magnitude and distance with radius. This 
suggests that a brightening bias due blending or crowding, 
rather than metallicity, could be the driving factor of the 
trend between radial distance and magnitude trends in M31. 


Figure 15. The distance determinations from Table[^are plotted 
here, separated by filter. In the left half of the plot, the DIRECT 
distances and the PHAT distances using the DIRECT Cepheids 
are presented for both Galactic and LMC anchors. The the right 
half of the plot, the PHAT distances in the B and I filters are 
shown for both Galactic and LMC anchors, as well as the NIR 
distance determination. 


7 DISCUSSION 

We present our distance estimates for the entire dataset 
(from Table in Table This dataset has the advantage 
of having the largest sample size of Cepheids in M31 in the 
visual bands from HST, as well as having a larger sample 
in the NIR bands than previous studies (Riess et al. 2012; 
Kodric et al. 2015). 

Table and Figure show that there are significant 
biases in the slopes and intercepts of the PL relations of 
the ground-based DIRECT survey compared to the PHAT 
survey. The observed differences suggest that Cepheid ob¬ 
servations obtained from ground-based surveys, while com¬ 
prehensive in scope, are not ideal tools in the era of precision 
cosmology, and that Cepheid observations of more distant 
galaxies with HST may suffer from similar effects. Surpris¬ 
ingly, the biases do appear to cancel somewhat, leading to 
only very modest changes in the mean distance modulus. 

In Table we give recent values from other studies 
for the distance modulus of M31. For a full compilation 
and discussion of M31 distances, we refer the reader to de 
Grijs & Bono (2014). The suggested M31 distance modu¬ 
lus from their study, incorporating various distance mea¬ 
surement techniques and statistically weighted errors, is 
24.46T0.10 (assuming an LMC distance modulus of 18.50). 
Our distance modulus of 24.51 ± 0.08 determined from the 
NIR PL relations is in agreement with this range. 

In contrast to the NIR results, the visual Wesenheit 
(BI) distance modulus of 24.32 ± 0.09 (an average of the 
MWG and LMC distances) falls outside of the 1-a range of 
de Grijs & Bono (2014). The distance moduli determined 
from the Ws/ relations for the complete Cepheid sample 
and Wvj relations for the DIRECT Cepheids with PHAT 
magnitudes are shorter than published distances, for both 
LMC and MWG anchors (see Table [^. They also disagree 
with our NIR distance determination of 24.51 ± 0.08. Al- 
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Table 8. M31 Distance Moduli 


Distance modulus 

Filter(s) 

Source 

Calibration 

Anchor 

Source 

24.44±0.10 

BVRI 

Ground-based 

Gepheids 

LMG 

Madore & Freedman (1990) 

24.44±0.12 

B, V 

Ground-based 

Eclipsing Binary 

- 

Ribas et al. (2005) 

24.32±0.12 

B, V 

Ground-based 

Gepheids 

LMG 

Vilardell et al. (2007) 

24.47±0.12 

1 

Ground-based 

TRGB 

- 

Durrell et al. (2001) 

24.47±0.07 

I 

Ground-based 

TRGB 

- 

McConnachie et al. (2005) 

24.49±0.11 

R, I 

Ground-based 

Gepheids 

LMG 

Joshi et al. (2003) 

24.41±0.21 

R, I 

Ground-based 

Cepheids 

LMG 

Joshi et al. (2010) 

24.32 ± 0.09 

B, I 

HST 

Cepheids 

This study 


24.5±0.1 

F606W, F814W 

HST 

RR Lyrae 

Carretta et al. (2000) 

Brown et al. (2004) 

24.38±0.05 

V, I (F555W, F814W) 

HST 

Cepheids 

LMC/Key Project 

Freedman et al. (2001) 

24.38 ±0.064 

IR (FllOW, F160W) 

HST 

Cepheids 

NGC 4258 

Riess et al. (2012) 

24.51 ± 0.08 

IR (FllOW, F160W) 

HST 

Cepheids 

This study 


24.54T0.08 

H (F160W) 

HST 

Cepheids 

LMG 

Macri et al. (2001) 


though our mean visual distance determination is within 1 
sigma of several recent studies (Freedman et al. 2001, Ribas 
et al. 2005, Vilardell et al. 2007, Joshi et al. 2010, among 
others), it does fall on the short side of the range of recently 
published distances, as well as slightly outside the determi¬ 
nation of de Grijs & Bono (2014). 

Unfortunately, there are no other similar published HST 
studies in comparable visual filters for Cepheids in M31 for a 
comparison. As examined by Mochejska et al. (2000), blend¬ 
ing in the visual bands in ground-based observations may 
cause underestimates in Cepheid distances by 6 to 9%. Al¬ 
though ground-based studies are more affected by these ef¬ 
fects than HST, if further blending or crowding effects re¬ 
main it could partially account for distances that remain 
shorter than expected. Williams et al. (2014) show that there 
is a brightening bias for the PHAT photometry in M31, es¬ 
pecially for the inner region and fainter magnitudes. For 
a typical Cepheid, the infrared bands and F814W magni¬ 
tudes can be biased up to 0.1 magnitudes (Williams et al. 
2014). The F475W magnitudes are less affected, exhibiting a 
brightening of <0.05 magnitudes; colors also appear to have 
minimal bias at the level of a few hundredths of a magnitude. 
Magnitudes may be brightened up to 0.1 in the inner regions 
compared to ~0.02 at 25 kpc from the center of M31. These 
biases reflect the crowding and blending problems present. 

In a more direct comparison, we can compare our dis¬ 
tance result in the infrared to similar NIR Cepheid studies in 
M31 with HST. Riess et al. (2012) obtain a distance modu¬ 
lus of 24.38±0.064. When we account for the increase of the 
distance to NGC 4258 by Humphreys et al. (2013) from 7.4 
to 7.6 Mpc, this pushes the Riess et al. (2012) value to a 0.06 
farther distance. This leads to a difference between the Riess 
et al. (2012) value and our NIR distance of l-cr. The differ¬ 
ence may be explained by differing global fits, where our 
slightly steeper slope and fainter intercept pushes towards a 
further distance. Our distance estimate from the long period 
(greater than 10 days) Cepheids using the scale from Riess 
et al. (2009, 2011) is also consistent with the updated Galac¬ 
tic and LMC distance calibrations done by Bhardwaj et al. 
(2015) in the infrared. The distance determination by Ko- 
dric et al. (2015) gives a 0.068 closer distance to M31 relative 
to the Riess et al. (2012) value. This puts their distance ap¬ 
proximately 1-CT lower than Riess et al. 2012 and within 2-a 


of our distance estimate, when accounting for the difference 
in the updated maser distance. Kodric et al. (2015) obtain 
their Wesenheit magnitudes slightly differently, making a 
direct comparison incomplete. Additionally, they compare 
their PL ht to Riess et al. (2012) at log(P) = 1.2, whereas we 
re-solve the system of M31 and NGC 4258 Cepheids. The 
hnal results of each study appear to be sensitive to either 
the included Cepheid sample, the distinct method for deter¬ 
mining the distance, or both. 

Derived from the complete sample of Cepheids, our NIR 
distance estimate of 24.51 ± 0.08 and visual distance of 24.32 
± 0.09 disagree slightly beyond the l-cr level; it is unclear 
from where this disagreement stems. We do not find evi¬ 
dence for bias from random sampling and the optical bands 
face less crowding, as shown in Williams et al. 2014. It is also 
possible that blending from nearby companions has affected 
the visual filters more than the infrared filters, causing the 
visual distance estimate to be shorter. Previous studies sug¬ 
gest that the contamination due to blending may be greater 
in the visual than in the infrared, though more rigorous stud¬ 
ies are needed to determine the effects and extent of blending 
on observations (Mochejska et al. 2000, Gieren et al. 2008). 

We are inclined to trust the NIR distance estimate over 
the visual band estimate for several reasons. Though the vi¬ 
sual bands are not in the same native photometric system as 
their calibrating relations from the Milky Way and LMC and 
the transformations to B, V, and I can be tricky. Although 
the NIR Wesenheit relation is also in a different photomet¬ 
ric system than the original Riess et al. (2009) and Riess 
et al. (2011) papers, it is in the same system as Riess et al. 
(2012). Additionally, it is commonly thought that the longer 
wavelengths of the infrared reduce the effects of extinction, 
amplitude, and metallicity on Cepheid observations (Madore 
& Freedman 1991, Persson et al. 2004). Due to the absence 
of significant scatter in the Wesenheit BI relation, extinc¬ 
tion effects are probably not severe. However, the remaining 
effects could introduce greater variations and possible biases 
in the visual bands than the infrared. Therefore, we suggest 
that our NIR distance estimate is a more reliable distance 
determination. 
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8 CONCLUSIONS 

We have analyzed Cepheid variables located in the PHAT 
data from the PAndromeda, DIRECT, and POMME (as 
available in Riess et al. 2012) ground-based surveys of M31. 
Through analyzing the visual and NIR magnitudes obtained 
from HST random phase observations, we make the follow¬ 
ing conclusions: 

1. We use a sample of Cepheids to construct a Period- 
Luminosity relation using highly accurate HST magnitudes 
obtained by the PHAT survey. In particular, the visual We- 
senheit relation from the complete sample (Wbi) shows less 
scatter over previous studies in the visual bands, leading 
to smaller random uncertainties in the distance modulus of 
M31. 

2. The dispersions in the visual Wesenheit and the NIR 
Wesenheit relations are very similar. The dispersion in the 
NIR relations is comparable to that derived by Riess et al. 
(2012), but with a 2.5x larger Cepheid sample that further 
reduces the uncertainty in the M31 distance modulus. Al¬ 
though our dispersion is slightly larger than that of a similar 
study by Kodric et al. (2015), they use a different calculation 
of the Wesenheit magnitude. 

3. We obtain a value of the distance modulus for M31 
of 24.32 ± 0.09 from the visual filters and 24.51 ± 0.08 
in the NIR using ACS and WFC3 photometry. These values 
are both consistent with recently published distance moduli. 
However, they disagree with each other by slightly more than 
l-cr due to the distance estimate at visual wavelengths being 
on the shorter end of published distances. 

4. The PHAT survey provides highly accurate HST 
magnitudes of 175 Cepheids in M31. The superiority of the 
PHAT photometry from HST is clear in this dataset from 
the significantly smaller dispersions than that found using 
ground-based data, despite the fact that the Cepheid magni¬ 
tudes are random phase. This is likely due to the enormous 
improvement in photometric quality and resolution as com¬ 
pared to ground-based surveys. 

5. We find a statistically significant magnitude radial 
trend leading to a distance modulus-radius relationship, 
which does not appear to be explained by a metallicity cor¬ 
rection but may be due to crowding or blending. Further 
work must be done to conhrm or rule out metallicity as a 
cause of distance discrepancies between inner field and outer 
field Cepheids. 
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